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Here, we report the draft genome sequence of the Mycobacterium tuberculosis-lilie Mycobacterium smegmatis mutant strain, 
mc^51, compared to that of wild-type strain M. smegmatis mc^l55. The draft genome sequence comprises 6,823,739 bp, revealing 
6,569 coding DNA sequences (CDSs) and 50 RNA genes (4 rRNA genes and 46 tRNA genes). 
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It is known that a successfully established infection of Mycobac- 
terium tuberculosis is partially due to its ability to survive and 
persist in macrophages (1). To defend against infection, the pro- 
duction of reactive oxygen species (ROS) by the host is one of 
important classical innate mechanisms. Accordingly, to overcome 
ROS stresses, mycobacteria have evolved many detoxification 
strategies to scavenge H2O2, the major contributor of ROS (2). As 
a model bacterium, Mycobacterium smegmatis is widely used for 
investigating mycobacterial physiology, drug selection, and vac- 
cine development (3-8). Historically, continuous selection was an 
effective method for making mutant strains. For example, Myco- 
bacterium bovis BCG, the only vaccine against tuberculosis, was 
attenuated by serial passage of M. bovis in vitro over 13 years (9). 
We used a similar adaptive evolution strategy that mimicked the 
host environment of H2O2 stress by gradually increasing H2O2 
concentrations to select H202-resistant M. smegmatis mutants. 
After 54-day selection, we obtained the M. tuberculosis-\ike M. 
smegmatis mutant strain, mc^51, with an 80-fold higher MIC of 
H2O2 and a 100-fold lower MIC (0.1 mg/liter) of isoniazid com- 
pared to those of wild-type M. smegmatis mc^l55; this is similar to 
M. tuberculosis with its high resistance to H2O2 and conversely, 
susceptibility to isoniazid. In addition, the mutant strain mc^51 
shows slow growth and has a growth advantage in macrophages 
compared to wild-type mc^l55. 

We therefore sequenced the whole genome of the mutant 
strain to explore the genomic variations that are potentially in- 
volved into its phenotype. The genome of M. smemgatis mc^51 
was sequenced on an lUumina HiSeq 2000 platform. A total of 
4,123,831 paired-end reads with a length of 100 bp, corresponding 
to 121-fold coverage of the genome, were generated. Raw reads 
were first filtered using the DynamicTrim and LengthSort Perl 
scripts provided in the SolexaQA suite (10) and then assembled 
using SOAPde«ovo (http://soap.genomics.org.cn), yielding 96 
scaffolds with an average length of 71,080 bp. The draft genome 
comprises 6,823,739 bp, with 67.4% G+C content and an N^g 
length of 150,887 bp. The genome was annotated using the Rapid 
Annotations using Subsystems Technology (RAST) program (11), 



revealing 6,569 coding DNA sequences (CDSs) and 50 RNA genes 
(4 rRNA genes and 46 tRNA genes). 

We then used the Burrows-Wheeler Aligner (BWA)/SAMtools 
(12) to call single nucleotide polymorphisms (SNPs). Compared 
with the parental strain mc^l55, 22 SNPs were identified in the 
mutant genome. Nineteen of these mutations are located in the 
coding regions, and 12 of the 19 SNPs are nonsynonymous muta- 
tions. Interestingly, we found nonsynonymous mutations in sev- 
eral genes encoding detoxified enzymes, including dehydroge- 
nase, oxidase, and monooxygenase. We also identified a mutation 
in the fur gene encoding a ferric uptake regulator, which is located 
immediately upstream of the catalase-peroxidase KatG. The SNPs 
of fur were also found in isoniazid (INH)-resistant clinical iso- 
lates, which implies that the mutants of Fur might be produced in 
humans. 

The information will undoubtedly shed light on the ways in 
which mycobacteria adjust their own fitness and acquire persis- 
tent in macrophages, which might give new insights into design- 
ing control strategies against M. tuberculosis. 

Nucleotide sequence accession number. The draft genome se- 
quence has been deposited at DDBI/EMBL/GenBank under the 
accession no. JAIDOOOOOOOO. The version described in this paper 
is the first version. 

ACKNOWLEDGMENTS 

This work was supported by grants from the National Natural Science 
Foundation of China (no. 31270178 and 31070118). 

REFERENCES 

1. Podinovskaia M, Lee W, Caldwell S, Russell DG. 2013. Infection of 
macrophages with Mycobacterium tuberculosis induces global modifica- 
tions to phagosomal function. Cell. Microbiol. 15:843-859. http://dx.doi 
.org/10.1111/cmi.l2092. 

2. Pieters ]. 2008. Mycobacterium tuberculosis and the macrophage: main- 
taining a balance. CeU Host Microbe 3:399-407. http://dx.doi.org/10.10 
16/j.chom.2008.05.006. 

3. Zahrt TC, Song 1, Siple J, Deretic V. 2001. Mycobacterial FurA is a 
negative regulator of catalase-peroxidase gene katG. Mol. Microbiol. 39: 
1174-1185. http://dx.doi.Org/10.llll/j.1365-2958.2001.02321.x. 



January/February 2014 Volume 2 Issue 1 e00092-14 



Genome Announcements 



genomea.asm.org 1 



Li et al. 



4. O'Toole R, Smeulders MJ, Blokpoel MC, Kay EJ, Lougheed K, Williams 

HD. 2003. A two-component regulator of universal stress protein expres- 
sion and adaptation to oxygen starvation in Mycobacterium smegmatis. ]. 
Bacteriol. 185:1543-1554. littp://dx.doi.org/10.1128/JB.185.5.1543-1554. 
2003. 

5. Venkatesh J, Kumar P, BCrishna PS, Manjunath R, Varshney U. 2003. 
Importance of uracil DNA glycosylase in Pseudomonas aeruginosa and Myco- 
bacterium smegmatis, G+C-rich bacteria, in mutation prevention, tolerance 
to acidified nitrite, and endurance in mouse macrophages. J. Biol. Chem. 
278:24350-24358. http://dx.doi.org/10.1074/jbc.M302121200. 

6. Tao J, Han J, Wu H, Hu X, Deng J, Fleming J, Maxwell A, Bi L, Mi 
K. 2013. Mycobacterium fluoroquinolone resistance protein B, a novel 
small GTPase, is involved in the regulation of DNA gyrase and drug 
resistance. Nucleic Acids Res. 41:2370-2381. http;//dx.doi.org/10.109 
3/nar/gksl351. 

7. Osborne R. 2013. First novel anti-tuberculosis drug in 40 years. Nat. 
Biotechnol. 31:89-91. http://dx.doi.org/10.1038/nbt0213-89. 

8. Sweeney BCA, Dao DN, Goldberg MF, Hsu T, Venkataswamy MM, 
Henao-Tamayo M, Ordway D, Sellers RS, Jain P, Chen B, Chen M, Kim 



J, Lukose R, Chan J, Orme IM, Porcelli SA, Jacobs WR, Jr. 2011. A 
recombinant Mycobacterium smegmatis induces potent bactericidal im- 
munity against Mycobacterium tuberculosis. Nat. Med. 17:1261-1268. 
http://dx.doi.org/10.1038/nm.2420. 
9. Calmette A, Guerin C, Boquet A, Negre L. 1927. La vaccination preven- 
tive centre la tuberculose par le "BCG." Masson et cie, Paris, France. 

10. Cox MP, Peterson DA, Biggs PJ. 2010. SolexaQA: at-a-glance quality 
assessment of lUumina second-generation sequencing data. BMC Bioin- 
formatics 11:485. http://dx.doi.org/10.1186/1471-2105-ll-485. 

11. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, 
Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson 
R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, 
Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, 
Wilke A, Zagnitko O. 2008. The RAST server: Rapid Annotations using 
Subsystems Technology. BMC Genomics 9:75. http://dx.doi.org/10.1186 
/1471-2164-9-75. 

12. Li H, Durbin R. 2009. Fast and accurate short read alignment with 
Burrows-Wheeler transform. Bioinformatics 25:1754-1760. http://dx.doi 
.org/ 1 0. 1 093/bioinformatics/btp324. 



2 genomea.asm.org 



Genome Announcements 



January/February 2014 Volume 2 Issue 1 e00092-14 



